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ABSTRACT 

Calculations for inclusive t/j production in T decay through the color-octet mechanisms 
ccp^+^Lj, 8] + g are presented. It is shown that these 0{a^v^) contributions 
compete with other color-octet mechanisms considered in the literature. A critical numer- 
ical analysis of the color-octet contributions to T — > -i/^ -f X shows that further work in this 
channel, both theoretical and experimental, is necessary in order to clearly understand the 
significance of the color-octet component of the cc component inside the system. 
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I. INTRODUCTION 

The unexpected large rates for if: and •0' prompt production at the Fermilab Teva- 
tron[l] gave rise to two interesting theoretical developments. On the one hand it was 
realized that, at high energies, the most important contributions to t/^ production comes 
from processes with a fragmentation interpretation [2]. On the other hand it was necessary 
to develop new ideas for heavy quarkonium summarized in an effective theory called non- 
relativistic QCD (NRQCD) [3]. This theory describes the interactions of non-relativistic 
quarks and is typically applied to QQ bound states such as upsilon or -0. The NRQCD 
Lagrangian is made precisely equivalent to full QCD through the addition of local interac- 
tions that systematically incorporate relativistic corrections through any given order in the 
heavy quark velocity v. The theory has a simultaneous expansion in Ug and in the velocity 
of the heavy quark. The size of a term in the NRQCD Lagrangian can be estimated using 
well defined velocity counting rules [4]. 

Creation and annihilation of heavy quarkonium are naturally described within this 
theory which incorporates the factorization ideas. In fact, these processes, which necessar- 
ily occur at very small distances fa ^ [mq denotes the heavy quark mass) as compared 
with the size of the meson ?a — —. are described by the short-distance coefficients entering 
in the relativistic corrections to the leading order NRQCD Lagrangian. The evolution of 
the quark pair into physical quarkonium involves nonperturbative effects which are en- 
coded within matrix elements of higher order operators with a well defined hierarchy in 
the V expansion. 

As a consequence of NRQCD, the possibility arise for the formation of physical quarko- 
nium through higher Fock states. In particular, we have the possibility for the formation 
of physical quarkonium from a QQ pair in a color-octet configuration. This is the most 
attractive explanation to the surplus V' and V'' production at the Tevatron [5]. 

Unfortunately the NRQCD matrix elements can only be computed from lattice calcu- 
lations which are still in their infancy. The way people proceed is to fit these parameters to 
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some experimental data and use their values to make predictions for other processes where 
they are involved. Using the values for the NRQCD matrix elements as extracted from 
the fits to charmonium hadroproduction [6,17], predictions have been made for color-octet 
signals in prompt quarkonium production in e^e~ colliders [7], Zq decays [8], photoproduc- 
tion at fixed target experiments [9], hadroproduction experiments[10] and B-meson decays 
[11]. 

Inclusive if) production in T decay is a very good place to test the NRQCD ideas 
[12]. Unlike other •0 production processes , the short-distance coefficients highly suppress 
color singlet contributions to T — > t/; -|- X. In fact, color-singlet contributions start at 
through the 66[^5i,l] cc\^Si,l\ + gg and 1] cc[3 ;Si,l] + gggg mecha- 

nisms whereas color-octet contributions to the same process start at 0{a^v^) through the 
66[^<S'i, 1] cc[^ 51,8] +gg mechanism. Hence, the short distance factors of the color-octet 
contributions are enhanced by factors of l/a1 as compared with the color-singlet ones. 
This enhancement easily overcomes the suppression of the long distance factors. 

In principle we also have 0{a^vf) contributions from the bb in a color-octet configu- 
ration. However, such contributions are naively suppressed by {vb/vc)"^ as compared with 
contributions from a cc in a color-octet configuration. Actual calculations refiect this naive 
counting suppression [12]. 

In reference[12] (^(ct^ff ) and leading electromagnetic contributions to T ^ i^+X were 
calculated within NRQCD. The most important mechanism turned out to be: 66["^5'i, 1] 
cc[^S'i,8] -I- gg^ which gives a BR ~ 2.5 x 10~^ when the value Q!s(2mc) = 0.253 is used 
in the numerical calculations as suggested by the fragmentation limit for the T — > V' + 
decay width normalized to the T ^ ggg decay width. The next largest contribution is 
the indirect production through XcJ decay considered by Trottier[13]. When the values 
for matrix elements extracted from the Tevatron data[6] are used, this contribution gives 
'^jBR{T —>■ XcJ + X —>■ i/j + X) ~ 5.7 X 10~^. There also exist a comparable indirect 
contribution from ip' having a branching ratio ~ 5 x 10~^[12]. Adding up all these and 
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other even smaller contributions a BR ~ 4 x 10~^ is obtained [12]. This result is within 2a 
of the CLEO data: BR^^p = 1-1 ± 0.4 x 10~^ [14] which suggest other contributions need 
to be computed in order to bring theoretical predictions closer to experimental data. 

Below I compute the 0{a^gV^) one-loop contributions to inclusive t/^ production in T 
decay through the color-octet mechanisms hh^Si,V\ cc[^^+^Lj, 8] + g. The importance 
of these mechanisms can be easily understood from the calculation for the electromagnetic 
contribution to the same processes[12]. This is an 0{al^asV'^) contribution and gives a 
BR ~ 1.6 X 10~^. An 0{a'^) contribution should give roughly a BR ~ 10~^ or even 
larger if we consider the e^e^ ~ 5 x 10~^ factor coming from the quark electromagnetic 
couplings. The color-octet processes studied in this paper are also interesting due to the 
distinctive signature that the has a sharply peaked energy distribution described by a 
delta function with a peak at 5.2 GeV and is recoiled by a gluon jet. The delta function 
must be understood in the sense of NRQCD calculations, i. e., as a narrow distribution 
whose width is of order rricv'^ in the ip rest frame. When boosted to the T rest frame the 
width decreases by a factor ~ 0.3, therefore giving narrow distribution peaked at 5.2 
GeV whose width is ~ 150 MeV. 

II. THE bb[^Si] cc[^'+^Lj,8] + g MECHANISMS 

The factorization formalism developed in Ref. [3] can be easily generalized to the case 
of inclusive charmonium production from bottomonium decay[13]. For the process under 
consideration in this work we have the factorization formula 

dr{r^iP + X) = J](ir(m,n)(T|0(m)|T>(0^(n)>, (1) 

m,n 

where dT{m, n) is the short-distance factor for a bb pair in the state m to decay into a cc 
in the state n plus anything. The subscript m, n denote collectively the color and angular 
momentum quantum numbers of the heavy quark pairs. Contributions that are sensitive 
to the quarkonium scales (Bohr radius or larger) and to Aqcd can be absorbed into the 
NRQCD matrix elements (T|0(m)|T), (O'^(n)) . The relative importance of the various 
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terms in the double factorization formula (1) depend on the order of v in the NRQCD 
matrix elements and on the order of ag in the short-distance factors. 

The Feynman diagrams for the 0{agV^) contributions to the inclusive production 
in upsilon decay coming from the color-octet mechanisms bb[^Si,l] cc[^^'^^Lj,8\ +g are 
shown in Fig. 1. Calculations for these processes are very similar to the calculations for 
the radiative decays of quarkonium: QQ[^Si,l\ gg[^*+^Lj,l] -|- 7 studied in reference 
[15]. 

Using the standard techniques for heavy quarkonium calculations [6,16,21] the follow- 
ing amplitude is obtained for a bb pair annihilating into three off-shell gluons: 

M[bWSi] ^ gi{a)g2{b)g3{c)] = TabcA, (2) 

where Tabc = is the color factor 

ki ■ {k2 + k^)k2 ■ {ks + ki)k3 ■ {ki + ^2) 

LI stand for the Lorentz invariant structure 

LI =ei ■ e2{ki ■ kse^ ■ k2e ■ ki + k2 ■ k^es ■ kie ■ k2 + ki ■ k^k2 ■ kse ■ £3) 

-I- e • €3(^1 • k2ei ■ kze2 ■ ks — ki ■ ksei ■ k2e2 ■ ka — k2 ■ kse2 ■ kiei ■ k^) (3) 
+ (ei, ki ^ €3, k^) + (€2, k2 ^ €3, A;3); 

here, e denotes the bb^^Si] polarization vector and (/cj) stand for the polarization vector 
(momentum) of gi. 

The amplitudes quoted in this paper are obtained by projecting the free quarks am- 
plitudes over definite angular momentum ( '^^'^^Lj) and color (i, 8) quantum numbers in 
the usual way [6,16]. 

The two gluon fusion into a cc[^*"'"-^Lj, 8] processes are described by 

M[gi{a)g2{b) 5o,8c] = ^^.^ ^£(ei, fci, €2, ^2), (4) 
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riahc 4 1 

M[g^{a)g2{h) Pj,8c] = :r^9!^{i-^) Aj, (5) 



where 

^0 = [{ki ■ k2 + 4m^)(ei • €2^:1 ■ k2 - ei ■ k2e2 ■ ki) + kjei ■ k2e2 ■ k2 

+ k^ei ■ kie2 ■ ki - kfklei ■ 62 - ei ■ kie2 ■ k2ki ■ /C2] 
Ai = mc[kle{e* ,€1,62, /C2) + ei • kie{e*,e2, ki, /C2) + {ki, ei} ^ {^2,62}] 
A2 = VSml [ki ■ k2ei^e2u + k2,^ki^ei ■ 62 - ki,^e2^ei ■ k2 - k2,^ei^e2 • kile*/^", 

here, e*(e*'*^) denote the polarization vector (tensor) for the outgoing spin-1 (2) bound 
state and e stands for the Levi-Civita tensor. 

Using the hehcity projector techniques described in reference [15], the helicity amph- 
tudes can be readily calculated for the different processes. A 1/2 factor is inserted in the 
amplitudes for the full processes in order to avoid double counting of the one-loop diagrams 
as the amplitudes in equations(2),(4)and (5) are symmetric in the two gluons involved in 
the loop. Performing the loop integration and matching the theories by the procedure 
described in reference [6] (alternatively the cross section for free quarks annihilation can 
be computed and matched with the NRQCD results by expanding in the quarks relative 
momentum and identifying the different terms by a procedure described in reference[21]) 
the following results are obtained for the short distance coefficients: 

T[bb['S,,l] ^ cc['So,8ig] = Fr{l - r)\H%r)\^ (6) 
r[bbfS„l] ^ ccfPj,8],9] = (2j+^i)^f (l-^)El^i'^(^)l' 

where 

5 5 af(T|Oi(3gi)|T) 
^= 32486 ' ^ = (-'=/-^)- 

The explicit formulae for the helicity functions H^"^"^ are rather length and are deferred to 
the appendix. 
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The decay width for T — > ip+X is obtained by adding up the short-distance coefficients 
hsted above multiphed by (T|Oi(^S'i)|T) and their respective (Og (^^"'"^Lj)) long-distance 
matrix elements: 

r(T cc[''+'Lj,8]g - ^ + X) = 4^^|;^<^l^^^'^^)l^>^ 

where 

The / factors appearing in equation(9) are obtained by evaluating Hj^"^^ air — {rric/mi,)'^ = 
0.118 : 

= r(l - r)\H'{r)\^ = 5.79, /(°) = 2r(l - r)H^'^\r) = 6.54, 

/(I) = 2r(l - r)H^^\r) = 7.32, J^^) = 2r(l - r)H^^\r) = 8.34, 

where 

H(^\r) = J2\Hi'\r)\\ (10) 

i 

Using the heavy quark symmetry relation (Og ("^Pj)) = (2 J + l)(Og('^Po)), equation (9) 
can be written in the compact form: 



B=Uf'(OtCSo)) + r^^^i^] (11) 



32L- \ "// w 

where 

2 

/P = ^ /(•^) = 22.20. 

For numerical analysis purposes I list the decay width for inclusive production in 
T decay through the main mechanism studied in reference [12] 



r(T ^ cc['S^,8]gg ^ + X) ^ 4§^^<^l^i^''^i)l^>(^8 ^'"^i)) ^ ^^'^^^^^ ^^2) 

C I) 
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3. NUMERICAL ANALYSIS 

The decay width for T — > i/^ + X in equations (8) and (12) depends on as and the 
NRQCD matrix elements (T|Oif <Si)|T> , (oJ(3Po)> , (0^(3<5i)> and (oti^So)). In 
particular it is very sensitive to the chosen values for ag , (Og (^Pq)) and (Og (^S'l)). 
The values for the color-octet matrix elements have been mainly extracted from hadropro- 
duction and photoproduction data [6,11,17,18], these values being largely affected by NLO 
and perhaps by higher- twist corrections [17,19,20]. For the color-octet matrix element 
(O^i^Si)) a fit to CDF data gives [6] (Oti^Si)) E [0.0045, 0.0087]Ge1/3 _ ^j^g 
sis of the same data including higher order QCD effects such as initial state radiation 
of gluons by the interacting partons[19] gives somewhat smaller values: (Og (^/S'l)) e 
[0.0033, 0.0046]Gey3. As to {0^{^Po)) and {0^{^So)), different combinations of these 
matrix elements have been extracted from different hadronic processes [6,11,17,18,19]. The 
values extracted for the same combination entering in different processes differ roughly by 
a factor of 2. In a recent analysis [20] which summarizes fixed target hadroproduction 
data[17], Tevatron data[6,19] and photoproduction data[18], the following range was ob- 
tained for the combination of these matrix elements entering in the inclusive ip production 
at Tevatron : (OtCSo)) +3{Ot{^Po))/ml e [0.01,0m]GeV^ . 

On the other hand, production in e'^e~ annihilation was recently used to extract 
the same combination of matrix elements [22]. This process is not sensitive to (Og (^^o)) 
but is highly sensitive to (Og (^Pq))- The values obtained for (Og ("^Pq)) are very stable 
under changes in the input values for as, charm quark mass and the color-singlet matrix 
element {Of{^Si)) entering in the fit, thus allowing us to strongly constrain the values for 
this matrix element: {0^{^Po))/ml e [0.72,0.76] x IQ-^GeV^. 

In summary the only color-octet matrix element which seems to be firmly established 
is (Og (^Pq)) and Eqs.(8,12) should be analyzed as a function of the remaining color-octet 
matrix elements. 

The values a,(2mc) = 0.253, (©^(^^i)) = O.OUGeV^, {r\Oi{^Si)\r) = 2.3GeV^ 
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were used in [12] in the numerical evaluation of equation (12). The value for(T|Oi(^S'i)|T) 
is extracted from the Upsilon leptonic width and is very similar to potential model calcula- 
tions, thus reflecting the suitability of the the non-relativistic description for bottomonium 
system. The value for as is suggested by the fragmentation limit of equation (12) nor- 
malized to r(T — > ggg) which gives three times the fragmentation probability for a gluon 
fragmenting into a i/j, Pg^^tp [12,5]. However, fragmentation is not a good approximation 
for the process under consideration and the value used for (Og (^5'i)) seems to be too large 
on the light of updated data[6,17,20]. 

A direct evaluation of equation (12) was not performed in reference [12]. Instead, 
authors evaluated this equation normalized to r(T — > ggg) and multiplied by BR{T — > 
ggg) which was assumed to be ~ BR{T — > light hadrons) = 0.92. This procedure reduces 
the direct evaluation of Eq.(12) by a factor of j-^, which is compensated by the large 
value used for (Og (^^i)). A direct evaluation of Eq.(12) using ag{2mc) and the central 
value (oti^Si)) ~ 0.0066GeV^ obtained in reference [6] gives a branching ratio ~ 2.6 x 
10"'^. This BR is reduced by a factor ~ | if the central value quoted in reference[20] for 
(Og (^ySi)) is used in the numerical evaluations. 

Using the values ^^^"^^ = 0.0074^6^3 and (OtCSo)) = QMlGeV^ as extracted 
from i/j production in e+e" annihilation [22] and as{2mc) — 0.253 the contributions cal- 
culated in this work gives a branching ratio BR ~ 2.1 x 10~^. The pseudoscalar process 
Eq.(6) accounts for fa 25% of the contributions in equation (8), thus the branching ratio is 
less sensitive to (Og (^S'o)) than to (Og (^Pq)) • Adding up the contributions from equa- 
tions(8) and (12) with other even smaller contributions ( indirect psi production through 
Xcj and ip' etc) calculated in Ref. [12] a total branching ratio BR ~ 6.2 x 10~^ is obtained 
when as{2mc) is used in the numerical calculations. However, this BR is very sensitive to 
the chosen value for 0: = . 
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4. CONCLUSIONS 

Summarizing, inclusive i/j production in T decay through the color-octet mechanisms 
^^[^'S'l, 1] cc[^^~^^Lj,8\ + g is considered. A critical numerical analysis of these 0{alv'^) 
contributions and other in the current literature is performed. The total branching ratio 
for T — > '(/' + X and the size of the calculated contributions strongly depend on the values 
chosen for ag and the color-octet matrix elements involved in the different mechanisms. 
Using as{2mc) as suggested by the fragmentation interpretation of the 0{asV'^) term [12] 



extracted from e'^e" annihilation[22] , the contributions considered in this paper give a 
BR ~ 2.1 X 10~^. Adding up this branching ratio with the contributions calculated in 
reference [12] a total branching ratio BR ~ 6.2 x 10""^ is obtained in good agreement with 
the CLEO data BR^xp = 1.1 ± 4 x lO'^ and consistent with the ARGUS upper bound 



In the light of the available experimental information, the calculations presented in 
this paper indicates that color-octet mechanisms account for most of Psi production in 
Upsilon decay. However, further work, both theoretical and experimental, is necessary in 
order to have a clear idea on the role of the color-octet component of the cc pair inside the 
ip system. 

From the experimental point of view, is necessary to remove the existing inconsisten- 
cies as the ARGUS collaboration results are not confirmed by the CLEO collaboration 
[14]. A measurement of the decay width and of the energy spectrum of the ip would be 
desirable as the later can discriminate between the different ip production mechanisms. In 
particular, ijj^s produced through the color-octet mechanisms considered in this paper have 
the distinctive signature that the t/^'s are sharply peaked in energy and are recoiled by a 
gluon jet. This is in contrast with other mechanisms where a spread distribution in energy 
is expected. 

From the theoretical point of view, the size of the calculated contributions strongly 



and the central values 




recently 



BRexp < 0.68 X 10-3 [14]. 
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depend on the assumed value for ag. A numerical evaluation of equations (8, 12) using 
asi^rrib) decrease the calculated BR roughly by a factor of \. Hence, a calculation of the 
0{a^gvl) color-octet mechanisms hb^Si,X\ ~^ cc^^^^Lj,'S\ggg and the 0{a^g) color-singlet 
mechanisms hh^Si,V\ — > cc^Si,V\gg and 66pS'i,l] — > cc^Si^'^gggg is necessary in order 
draw definitive conclusions. The color singlet contributions have been "crudely estimated" 
to give a branching ratio of a few xlO~^ [13]. 
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APPENDIX 

H'(r) =- [{Di{2x) - Di{0) - ^ ln{2x) - ^-^{2Di{x) - 2Di{0) + hn^{l - x)) 

=\/| + (1°^ + 4i^in(2))to(l - .) - -i^^lnHl - .) 

+ K + 2 , A ln{2x) + 3- ^ L>z 3 D%{2x) 

^x^ a;(l — 2a;)^ a;'^(2 — a;) x^ 

2 ^X X _^ , , 1 *^ 7 / \ "I 

— 4 jTT- :—D%[x) + iDTT— rm(l — x) . 

a;2(2-a;) ^ ^ a;(2 - a;) ^ 

^(i)(r) =A [l(i:)z(0) - i:»z(2a;)) - 6(1 - x){Di{x) - Di{2x) 

X^ X 

2 2 — 8a; + 

- Ini2)lnil - a;)) + -(1 - x)il - 2x)lnil - a;) + ^ Zn(2a;)l . 

X 1 — 2x 



H[^\r) := ^^\ ^ [- {Di(0) - Di(2x)) -x- 2x{Di(x) - Di{2x) - ln{2)ln{\ - x)) 

X X 

2 — X — 2a;^ 

+ Ml - a;) + 2(1 + a;)Zn(2a;)l . 

X 

Hi\r) J-4 [(6 - 5x)x + 2 ^^i?^±i?£!(i _ ,),„(i _ 

1 10-12a; + 5a;2 , .^^^ ^^ 2 6 - 38x + Tlx^ - 37a;^ , 

- (Dz(0) - Dz{2x)) + /n(2a;) 

d-a;")^ 1 

- 8 \ (Di(2x) - 2Di(x) - -In^il - x) + Di(0) + Mnil - x)) 

a;^(2 — a;) 2 

4 6 — 6a; — a;^ /, , , ttx 

+ 3 (M2)-«2) 

- ^ (12 - 26a; + ISa;^) {Di{x) - Di{2x) - ln{2)ln{l - x))] . 

iff )(r) JbHEl. \ _ 1(38 - 9a;)a; - -(4 - 13a; + 16a;^ - 4a;^)/n(l - a;) 
a;-^ 3 X 

+ 2^^:£1^ (Dz(0) - Di{2x)) - Y^i'^ - 11a; + 16a;2 - 4a;3)Zn(2a;) 

Zi X -L ZiX 

+ 8 (lufKi^^Jtf!) (i;^(2a;) - 2Di{x) - \ln^{l - x) + Di(0) + mln{l - x)) 

X [Z X j Zi 

- 4r (in(2) - 4) + 4(8 - 12a; + 3a;2) (Di(a;) - L>i(2a;) - ln{2)ln{\ - 

3 a; 2 
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^^2) (-^^ ^l^x+^{l-6x + 6x^)ln{l -x)- (5 - 6a; + 2x'^) {Di{0) - Di{2x)) 

x'^lz — x) 2 
+ - "-^^-^^-^ (/^(2) - - 16(1 - x){Di{x) - Di{2x) - ln{2)ln{l - x)) 

+ 4(1 -6a;)Zn(2a;)]. 
where x=l—r and Di denotes the Dilog function D%{s) = J* ^-^^^dt. 
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FIGURE CAPTION 

Fig.l. One of the six diagrams for the short-distance processes bb[^Si,l] — > ccf^^+^Lj, 8] + 
9- 
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